Raman thresholds and rigid to floppy transitions in calcium silicate glasses 
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Alkaline earth silicate glasses xCaO — (1 — x)Si02 exhibit a well marked threshold in Raman 
lineshapes which can be related to the onset of network rigidity as the concentration of calcium 
oxide x is decreased. The present results are analyzed by constraint counting algorithms and more 
deeply characterized by a size increasing cluster approximation that allows to perform Maxwell 
mechanical constraint counting beyond the usual mean-field treatment. This permits to discuss 
under which structural conditions an elastic intermediate phase can be obtained. 
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I. INTRODUCTION 

Magmatic liquids are the principal agents of mass and 
heat transfer in the Earth and terrestrial planets and 
intensive research has been accomplished to understand 
the processes of mass or energy transfer with respect to 
melt or structural properties p], • Viscosity, thermal or 
electrical conductivity are indeed directly related to the 
structure of silicate melts which furthermore control the 
temperature behavior of magmas and their formation or 
crystallization 

In this context, calcium silicates of chemical formula 
xCaO — (1 — x)Si02 have received little attention as com- 
pared to the corresponding alkaline systems so that much 
of their properties over the whole calcium glass formation 
range are still object of active research. Most studies 
have indeed focused on the x = 0.50 molar concentration 
which corresponds to the crystalline wollastonite compo- 
sition p. Several studies have stressed the special role 
played by the calcium atom which acts as a modifier in 
the silicate networks |6j and leads to a global increase of 
the density from 2.75 g.cm~ 3 at x — 0.38 to 2.97 g.cm~ 3 
at the x — 0.60 pyrosilicate composition 0,0- Extensive 
studies have been performed to understand the misci- 
bility limits in this system and in other alkaline earth 
silicates |9L fTo| . 

A special attention has been devoted to the coordi- 
nation number of the calcium atom using EXAFS 
and X-ray spectroscopy 0] . More recently, Yannopoulos 
and co-workers have been studying inelastic light scatter- 
ing of calcium silicates around the x = 0.50 composition 
by polarized Raman spectroscopy |13|. The results with 
composition show a marked change in behavior around 
47% calcium. Specifically, intensity ratios of particular 
Raman lines as well as the Boson peak frequency present 
a sharp jump at this concentration. On the other hand, 



a line at 606 cm~ x shows an abrupt increase but there 
remain some uncertainetv about its attribution to the 
so-called D2 ring line However, the conclusion ob- 
viously suggests the presence of a transition that has not 
been characterized by the authors. 

In this work, we show from the investigation of Ra- 
man spectroscopy that a very particular elastic state 
is reached in the glass when the concentration of cal- 
cium x equals 47%. This threshold is identified with a 
rigid to floppy transition from Maxwell mechanical con- 
straint counting and suggests that elastic transitions can 
take place in calcium silicate glasses as in chalcogenides 
|15l llq . In the latter, the addition of cross- linking units 
such as germanium or silicon into a basic network made of 
two- folded chalcogenide atoms (sulphur, selenium) con- 
straints indeed internal degrees of freedom of the network 
by increasing the number of bond-bending and bond- 
stretching forces that can lead to a very peculiar situa- 
tion when the number of constraints per atom equals the 
number of degrees of freedom • It has been identified 
by M.F. Thorpe from numerical simulations on amor- 
phous silicon as being a floppy to rigid transition which is 
characterized by the vanishing of the number of normal 
(floppy) modes of the dynamical matrix 0, 0|. The 
glass transition of such systems undergoing a rigid to 
floppy transition is also substantially affected |20|. Un- 
derstanding to what extent the methods used with suc- 
cess for the chalcogenides can be applied to oxide glasses 
is therefore not only of fundamental interest but also an 
attractive perspective for applied purposes jU ■ 

We first display the results of the Raman analysis and 
the deconvolution of the spectra which show a marked 
change in behaviour for some modes at x — 0.47. Next, 
we apply on this system mean-field Maxwell mechanical 
constraint counting that permits to compute the concen- 
tration x c at which the fraction of zero frequency (floppy) 
modes vanishes [HI ] . We improve the approach in Section 



2 



TABLE I: Glass transition temperatures of the (1 — x)SiC>2 — 
xCaO system. 



Concentration x 


42 


44 


47 


50 


53 


T 3 [°C] 


769 


769 


770 


781 


795 



IV by analyzing the present system with size increasing 
cluster approximations (SICA) that infer the effect of 
medium-range order on the nature and the location of 
the transition which is found to be in harmony with ex- 
perimental findings. Finally, we show which conditions 
in terms of medium range order have to be satisfied in 
order to obtain an intermediate self-organized phase 1221 
that is bounded by a rigidity and a stress transition [22| , 
in close correspondance with the chalcogen analog. 



II. EXPERIMENTAL RESULTS 
A. Sample preparation 

The samples were prepared by mixing pre-dried S1O2 
(99.99 %) and CaC0 3 (99.95 %) powders in the correct 
proportions. For each sample, the mixture was melted 
in a platinum crucible at 1650°C for two hours, and 
quenched by placing the bottom of the crucible in cold 
water. The samples were then annealed at a temper- 
ature around 760°C for five hours and cooled slowly to 
room temperature. The glasses were transparent and free 
of crystallization as confirmed by the absence of Bragg 
peaks in XRD spectra. Chemical microanalysis using 
the energy-dispersive X-ray (EDX) performed on opti- 
cally polished glasses showed a very low departure (less 
than 0.4 mol of Ca) between the theoretical and the real 
composition for all the studied samples. The glass tran- 
sition temperatures were determined with a differential 
scanning calorimeter Setaram DSC-1600 at a heating rate 
10°C /min. These values (see Table HJ are slightly larger 
than those reported by Shelby[24| using the dilatometric 
technique, but they exhibit the same global trend 



B. Raman scattering 

Raman spectra were obtained on a Jobin-Yvon T64000 
spectrometer, with CCD detection, and through a BX40 
Olympus microscope (objective xlOO). The excitation 
wavelength was the 514.532 nm argon line of a Coher- 
ent Innova 70 Spectrum laser. The power was typically 
300mW at the laser output (20mW on the sample). The 
low- frequency part range (< bOOcm^ 1 ) was obtained in 
triple substractive mode (gratings I800gr/mm). For the 
high frequency part (above 300 cm -1 ), the spectrometer 
was in simple monochromator configuration with Notch 
filter (grating 600gr/mm). This allows a sufficient re- 
covering spectral range for merging the spectra. All the 
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FIG. 1: Raman lineshapes for various calcium concentrations. 
The low wavenumber range was acquired in triple substractive 
mode (see text) which explains the higher noisy level. The 
spectrum of silica at the bottom of the figure is displayed for 
comparison. 



data on the different compositions were collected in ex- 
actly the same conditions, one following each other, in 
order to be the most confident possible in any variation 
between subsequent spectra. 



C. Peak deconvolution 

The experimental spectra were first corrected from 
first-order Bose-Einstein factor, and also from scatter- 
ing law (/ ~ v^). The latter (a second-order correction 
compared to Bose-Einstein factor) is not of primary im- 
portance but is theoretically needed to extract the Ra- 
man susceptibility, which is the relevant physical quan- 
tity which can be used to extract information on dy- 
namics. The Bose Einstein reduction does not include 
the prefactor 1 /to generally used for describing the low- 
frequency part of Raman spectra in glasses [25j , accord- 
ing to the Shuker-Gammon formalism |26| . Examination 
of the bare spectra shows that the main variations upon 
increasing Ca content from 0.42 to 0.53 lie at wavenum- 
bers above SOOcm^ 1 and then justify to focus attention 
on this spectral range, dominated by intratetrahedron 
or Si — O — Si motions. The boson peak reflects or- 
der at larger distances. After this, the spectra were re- 
constructed with the FOCUS software [23 , by using a 
log-normal law for the boson peak at low frequencies, 
and gaussian shapes for higher frequency ones. The low- 
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FIG. 2: Fits of the Raman spectra of the x = 0.44 sample. 
The numbers (1-9) correspond to the various gaussian decon- 
volutions of the spectra (see text for details). BP is the Boson 
peak and LB the luminscence background. 



frequency (boson peak) part was fitted to insure a cor- 
rect description of the mid-frequency range, due to its 
long tail, but is clearly out of the scope of the present 
paper. Moreover, recent hyperRaman measurements in 
silica tend to show that the Boson peak responsable for 
the excess density of states is directly observed in hyper- 
Raman and inelastic X-ray scatteing, meaning that the 
low-frequency Raman component implies other degrees of 
freedom and is then much more complicated to interpret 
|28j. For higher frequency modes, gaussian shapes gave 
a better reconstruction of the spectra than lorentzian or 
damped harmonic oscillator ones 
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III. DISCUSSION 



A. Spectra description 



FIG. 3: Compositional behaviour in wavenumber, linewidth 
and relative integrated intensity for the lines 4 and 5 of the 
Raman spectra in the (1 — x)SiC>2 — xCaO system. 



Figure^displays the experimental Raman data for the 
six Ca concentrations 0.42 to 0.53. Figure [5] exhibits the 
typical agreement between experimental and fitted spec- 
tra (Bosc-Einstein and scattering corrected). Results of 
the fitting procedure are displayed in the different panels 
of Figure |3 and 0] (line wavenumbers, linewidths and in- 
tegrated intensity). The whole concentration range was 
described with 9 modes (plus the boson peak) : 3 modes 
in the range below 550 cm -1 , 3 modes to describe the 
sharp and assymetric peak near 600 cm -1 , and 3 modes 
for the higher-frequency group around 1000 cm -1 . A 
broad luminescence background was needed to describe 
the intensity increase up to high relative frequencies, a 
feature that can be connected with the presence of 3d 
impurities in the Ca precursor. 

For the 1000 cmr 1 group, it is sometimes described 
by four components, or even more. We consider that, as 
only three bumps are unambiguously visible on the spec- 



tra, the lowest number of components giving satisfactory 
reconstruction of the experimental data is the best pos- 
sible choice. The four components needed by Kalampou- 
nis et al. look obviously necessary on their data for 
x > 0.55 and x < 0.40, but not in the range investigated 
here. 

Modes 2-3 are somewhat weak and large, and then are 
rather inaccurate (at least compared to other ones) as 
their possible parameter variation in the Ca range is hin- 
dered by experimental uncertainty. Mode 1 is discussed 
below. For the following modes 4 and 5, the accuracy is 
better. Both wavenumbers monotonously increase upon 
increasing Ca concentration (see Figure |3J). Their inten- 
sities show a transfer from 5 to 4 somewhat above 47% 
Ca. For linewidths, mode 4 smoothly broadens whereas 
mode 5 suddenly narrows above 50%. The greater ef- 
fects concern the higher frequency modes, 6, 7, 8 and 
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9. One can first consider the latter three modes, clearly 
connected in one broad band. This assignement slightly 
differs from 0], where mode 7 is assigned to Q° species, 
and from Frantz and Mysen [2{j where Q 3 corresponds 
to a weaker mode and mode 8 is assigned to Si-bonding 
oxygen stretching vibrations. The superscript n in Q n 
denotes the number of bridging oxygens (not connected 
to Ca) on a SiO^/2 tetrahedron. According to Zotov [30| . 
mode 9 can be assigned to Q 3 species, mode 8 to Q 2 
and mode 7 to Q , The wavenumbers for modes 7 and 

8 are practically constant. It is not the case of their in- 
tensities, which show a continuous increase for modes 7 
and 8, whereas mode 9 is rather constant up to 47% Ca, 
and decreases rapidly after. The wavenumber of mode 

9 displays a decrease up to around 47% Ca, and then a 
slight increase. This effect is small but clearly larger than 
experimental error. One can note that the sum of inten- 
sities of modes 8 and 9 is roughly constant above this 
value 47% Ca (Figure |3}. For linewidths, modes 7 and 8 
are slightly increasing, whereas mode 9 is constant up to 
47% Ca, and shortens after (slightly, but clearly above 
experimental error) . One can then conclude on this part 
that the increase of Ca acts in two ways : a first one, 
continuous, leading to increase the number of Q 2 and Q 1 
species, to equilibrate the higher number of cations, and 
a second one, sharper, which consists in a transfer from 
Q 3 species to Q 2 ones upon increasing Ca content above 
47%. Then the remaining Q 3 are more decoupled of the 
network, as shown by their lower linewidth. 

Mode 6 is the most affected by the modification of Ca 
content as shown on Fig. 0] Its frequency is hugely in- 
creased (more than 50 cm -1 , which is considerable), the 
width of the line is decreased by more than a factor 2, 
between 42 to 48% calcium, and its integrated intensity 
falls down by one order of magnitude in the same con- 
centration range. When comparing with Ref. 13] where 
the 600 cm~ 1 asymmetric component is described only 
by one line, the present work shows that the anomaly 
standing around x = 47% is only due to mode 6, the 
other modes 4 and 5 overlapping in this component evolve 
monotonously in this composition range. This mode 6 
cannot then be assigned to Q° species as the correspond- 
ing intensity must increase upon enriching Ca content. 
It is more plausibly connected with Si — O — Si motions. 
The intensity dependence which is more rapid than for 
mode 9, would lead to think that mode 6 could corre- 
spond to Si — O — Si vibrations of Q 3 — Q 3 species. 

Mode 1 displays two specific features around 46 — 
47% Ca. Its frequency exhibits a noticeable downshift 
(20 cm -1 ), the width decreases by 30%, and the inte- 
grated intensity falls down by a factor two (Figure [SJ- 
Even if this mode is broader than the preceding ones, 
these parameter evolutions are sufficiently large to be un- 
ambiguous. This frequency downshift and line narrow- 
ing upon increasing Ca content is different from high- 
frequency modes (6, 9). This softening and narrowing 
is an uncommon behavior and can be qualitatively ex- 
plained only by a significant change of eigenvectors. Due 
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46 48 50 
Ca concentration % 

FIG. 4: Compositional behaviour in wavenumber, linewidth 
and relative integrated intensity for different lines (6-9) of the 
Raman spectra in the (1 — x)Si02 — xCaO system. The error 
bars are of the size of the symbols except for mode 6 and 9 
wavenumbers. 



to its frequency (it is the lowest-frequency mode apart 
from the boson peak), it implies some Ca vibrations. 
These vibrations would then more decoupled, and softer, 
above ca. 46%, that can be attributed to floppy regions 
in the glass. 



B. Raman threshold and rigid to floppy transitions 

The results obtained from Raman scattering and the 
behaviour with concentration of some modes closely par- 
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FIG. 5: Compositional behaviour in relative wavenumber, 
linewidth and relative integrated intensity of the mode 1 line 
lying at 392 cm" 1 for 42% Ca. 



allels the one found in chalcogenides |3l| ■ The frequency 
increase exhibited by modes 6 and 9 above 47% calcium 
can indeed be discussed in the context of rigid to floppy 
transitions. In fact, these modes look to harden just when 
the network becomes softer. This apparent contradiction 
can be explained by a non-uniform distribution of stress : 
in the stressed rigid range close to our observed threshold 
at x = 0.47, the network is mainly hardened by intercon- 
nected Q 3 species. Upon undergoing a rigidity transition, 
these links should become much less numerous, leading to 
a network which is dominated by Q 2 connections. The re- 
maining Q 3 islands should become somewhat decoupled 
from this network, producing an increase of the mode 9 
wavenumber and a noticeable line narrowing : the dis- 
tribution of Q 3 force constants is less important, leading 
to a linewidth closer to a lorentzian one. In the floppy 
range at high calcium concentration, the structure should 
be a network of Q 2 elements, separated by small 'harder' 
islands. The narrowing exhibited by mode 5 above 50% 
can be also interpreted in the same way. 



IV. ANALYSIS FROM CONSTRAINT THEORY 

Maxwell constraints counting appears to be useful to 
understand the present results as the calcium silicate net- 
work can be described by a molecular system constrained 
by bond-stretching and bond-bending (angular) forces. 



A. Maxwell global constraint counting 

We consider the CaO — S1O2 system as a network of 
N atoms composed of n r atoms that are r— fold coor- 
dinated. Enumeration of mechanical constraints ^?|]~ C3 
associated with bond-stretching (radial) forces leads to 
n c = r /2j while the number of bond-bending (angular) 
constraints is n@ — (2r — 3). The average number of 
floppy modes per atom F/N in this three-dimensional 
network is given bv|l9|: 

f = F/N = n d -n c = 3--^> r ( y -3) (1) 

r>2 

where is the dimension of the network and n c = n" + 
. Applied to the system of interest leads to: 

11 — 7x 

f = F/N = 3-^—^- (2) 

The latter equation holds if one assumes that silicon is 
four- fold, calcium and oxygen are two- fold coordinated 
[3^| . As one can see, the number of floppy modes vanishes 
when the network attains the critical concentration x — 
x c = 0.50 which is in close agreement with the thresholds 
observed from our Raman results. 

Equation J3J) defines a mean-field transition in which 
the number of mechanical constraints is computed from 
the macroscopic concentration x. Obviously, this elas- 
tic transition may be attained at x < x c provided that 
some macroscopic floppy subregions can emerge with the 
addition of alkaline earth oxide. On the other hand, it 
has been recently shown [1^|,[33 that the underlying na- 
ture of the floppy to rigid transition was more subtle 
and could contain under certain circumstances two 
transitions instead of the single one predicted by mean- 
fied constraint counting. In this context, stress will not 
spread randomly over the whole network as initially be- 
lieved but will accumulate in underconstrained subre- 
gions leading to the occurence of an intermediate elas- 
tic phase that is found to be stress free. The detection 
of this intermediate phase has been mostly accomplished 
from calorimetric probes using temperature modulated 
differential scanning calorimetry (MDSC) [34[. Here, a 
sinusoidal variation is added on the usual linear DSC 
ramp and permits to decompose the total heat flow into 
a reversing part that tracks the initial modulation and a 
residue. The latter is found to vanish in the intermediate 
phase H3. 

B. Cluster construction 

Size increasing cluster approximations (SICA) appear 
to be a useful tool to describe the elastic nature of the 
network backbone (floppy, intermediate, stressed rigid) 
and permits to take into account medium-range order 
effects that may, or may not, serve as ingredient for 
the presence of an intermediate phase. Furthermore, 



SICA start from the mean-field description as basic level 
and takes into account non-random structural elements 
and their related mechanical constraints in a system- 
atic fashion. This method has been first introduced to 
study the formation of fullcrcnes |35| or Penrose tilings 
in quasicrystals j^] and was then applied to quantify 
the boroxol ring statistics in amorphous B2O3 [37]]. Re- 
cently, SICA has been used in the context of floppy to 
rigid transitions |33| in an archetypal chalcogenide net- 
work (Ge x Sei- x ) and has led to the definition of a stress- 
free intermediate phase that depends substantially on the 
fraction of small rings in the structure. 

The basic level (/ = I ) of the SICA construction corre- 
sponds to the mean-field approximation, having as el- 
ements structural species that depend directly on the 
macroscopic concentration. The construction permits to 
generate clusters at step I = 2 sharing all possible com- 
bination of the basic (Z = 1) elements, clusters at step 
Z = 3 having three (Z = I) units, etc. The probabilities 
of the clusters are computed within the Canonical En- 
semble having energy levels E n related to bond creation 
between the basic level molecules. The construction is 
furthermore supposed to be performed at the formation 
of the network when T equals the fictive temperature Tf 
which is defined by the intersection of the extrapolated 
supercooled liquid and glass curves [38| . Mathematically, 
these probabilities will involve statistical weights (or de- 
generacies g(E n ) of a corresponding energy level) which 
correspond to the number of equivalent ways to connect 
two (Z = 1) basic units, and a Boltzmann factors [3^ of 
the form e„ = exp[— E n /TA. 

We have used as basic units the CaO molecule and 
a SW4/2 tetrahedron, having for respective probabilities 
the macroscopic concentration x and (1— x). The energy 
levels are defined from the consideration of all possible 
connections at the step Z = 2 and permit to distinguish 
the mechanical nature of the underlying clusters (floppy, 
isostatic, stressed). The creation of a chain-like Ca^Oi 
floppy structure [n c < 3) is related to an energy gain 
of Ef while the creation of a CaO connected to a sili- 
con tetrahadron (an isostatically rigid CaSiOs cluster, 
n c = 3) is associated with an energy Bj. Finally, the ba- 
sic network former represented by two connected SW4/2 
tetrahedra (a stressed rigid S^C^ cluster, n c > 3) cor- 
responds to an energy gain of E s . At step 1 = 2, three 
different clusters can be obtained (see fig. 0) and their 
probabilities are given below. We mention that edge- 
sharing SiO^/2 have been excluded from the construction 
as they have no experimental evidence at all in silicates 
[39|. The probabilities at setp Z = 2 are: 
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FIG. 6: Fhe different clusters obtained at SICA step I = 2. 
a) A fragment of the silica network with g(E s ) = 16, b) a 
wollastonite-like cluster CaSiOs with g(Ei) = 16 and c) a 
calcium-rich cluster (c) with g(Ef) = 4. 



P3 



4x 2 ef 



Ax2e f + 16x(l - x)e l + 16(1 



(5) 



out of which can be computed the (Z 
of calcium atoms 



,(2) 



P2 + 2pi 
4 - p 2 - 2pi 



2) concentration 



(6) 



Due to the initial choice of the basic units, the energy Ei 
will mostly determine the probability of isostatic clusters 
since the related Boltzmann factor e, is involved in the 
probability 10} of creating the isostatic CaSiO^ cluster (a 
CaO — bonding). In the case where E^ <C Ef,E s , 

the bonding of the network construction will be mainly 
achieved by isostatic clusters. 

For larger steps (Z > 2), one has to take care of possible 
isomers produced by two distinct clusters at the lower 
level. The cluster displayed in Fig [7^ can for example 
be produced either by connecting two tetrahedra 
onto a CaSiOs (Z = 2) cluster or by adding CaO and 
SiO^/2 onto a S12OA (Z = 2) cluster. The statistical 
weight g(E n ) of such isomers will be of course larger than 
the one of low connected clusters having a single pathway 
of production. 

As seen from equs. the cluster probabilities 

will depend only on two energetical parameters (i.e. the 
factors e s /ei and e//ej) and the problem can still be 
reduced by using a charge conservation law f° r tne 
calcium atom applied on the population of clusters of size 
Z. 



(7) 
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The ^-dependence of either e s /ei or e//ej means that 
either the Active temperature Tf or the energies E n de- 
pend on x but here only the e„ (a;) / e\ (x) dependence 
is relevant for the analysis. 

In order to obtain some clusters having significant 
medium range order, the SIC A construction has been 
realized up to the step / = 4. 



FIG. 8: Q n distribution in (1 — x)SiC>2 — xCaO as a func- 
tion of the calcium concentration, computed at SICA step 
I = 4. The vertical broken lines show the boundaries of the 
intermediate phase (see below). 



V. RESULTS 



C. Maxwell cluster constraint counting 

Next, one can apply on the generated set of clusters 
Maxwell constraint counting by enumeration of bond- 
bending and bond-stretching constraints and calculation 
of the corresponding expressions of n" and . Of par- 
ticular importance are the structures containing a ring 
having less than six atoms (see Fig. UJi), because one has 
to remove some extra constraints 

For each step / has been computed the total number 
of constraints 



,(0. 



Z^=i n c (i)Pi 

E^'i N iPi 



(8) 



where Mi is the total number of clusters of size I and n c ^ 
and Ni are respectively the number of constraints and the 
number of atoms of the cluster of size I with probability 
Pi. Applied to the set of clusters at step 1 = 2, one 
obtains for example: 



(2) 



22pi + 15p 2 + 4p 3 
6pi + 5p2 + 2p 3 



(9) 



We have determined either e s /ei or e//e, as a function of 
concentration, by solving the charge conservation law J7J). 
With these factors depending on the concentration x, it is 
possible to compute the cluster probabilities pi of a given 
step I as a function of composition and finally obtain the 
composition x c where the number of floppy modes fi = 
3 — Uc vanishes. It can be also computed the statistics 
of Q n species with Ca concentration as discussed next. 



A. Structural properties and speciation 

In this section, we consider the solutions of the SICA 
construction in terms of structure. One principal objec- 
tive of the present investigation is first to compare the 
model results with some experimental data on calcium 
silicate glasses such as the relative abundances and mix- 
ing properties of the structural units with respect to the 
concentration. 

Figures 00 and [§] represent the distribution of Q" units 
computed from SICA at step I = 4, respectively as a 
function of Ca concentration and as a function of the 
relative abundance of non-bridging oxygens NBO/T that 
is defined by: 



J2 n kx(Q k ) = NBO/T 



(10) 



where is the fraction of non-bridging oxygens per 
tetrahedra and x(Q k ) is the concentration of Q k 
species. The quantity defined by equation 1)10(1 is com- 
monly used to quantify network depolymerization with 
non-bridging oxygens |4lj . 

As one can see, the probability of finding a Q A unit 
decreases smoothly with the addition of calcium whereas 
the emergence of S^Og - (Q 2 ) and Si0 2 ~ (Q 3 ) groupings 
is noticeable. One can furthermore remark the crossover 
of the abundance of Q 2 and Q 3 at the concentration 
x = 0.47, a behaviour that has been also observed from 
the integrated intensity of mode 8 and 9 in the Raman 
spectra. 

Several limitations of the approach can however be ob- 
jected. First, the fact that we limit the cluster construc- 
tion to the step I = 4 does not permit to generate Q° 
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FIG. 9: Q n distribution in (1 — x)Si02 —xCaO as a function 
of NBO/T, computed at SICA step I = 4. The insert shows 
the calculated equilibrium constant Kci (solid line) and Kc2 
(broken line) as a function of NBO/T. See text for details. 
Data points are taken from Frantz and Mysen [29] . 



as these species are made of a central tetrahedron that 
shares four Ca atoms (that structure would be created 
by I = 5). As a consequence, the distribution of our Q n 
(n — 4, 3) are slightly overestimated in the high calcium 
range. Secondly, as the energetical factors are introduced 
only on the basis of the overall mechanical character of 
the bonding type (stressed, isostatic, floppy) between two 
basic units, difficulties arise in order to describe the be- 
haviour with changing alkaline earth cation as the size 
(or energetics, or steric hindrance,...) does not appear 
in this approach. This means that other ingredients [3^ | 
are necessary to describe the differences arising in Mag- 
nesium or Barium silicates. 

Next, we can focus on equilibrium constant |4"3 | K c \ 
and K c2 between species respectively related to the equi- 
libria: 



Si 2 Or ^ SiQ 2 + SiOi 



4SiOl~ ^ Si 2 2 5 ~ 



Si 2 O e 7 ~ 



(11) 



(12) 



The mass action constants can be computed from the 
above equilibria using the computed distribution of Q n 
species. It is obvious that eau. Hll|) will be the dominant 
equation in the low calcium region while both equilibria 
will have to be taken into account in the concentration 
range lying around the rigid to floppy transition when 
NBO/T is larger than the value 1. As one can see from 
the insert of Fig. HO the calculated K c \ from the equilib- 
rium constant of is in fair agreement with the two 
data points number reported by Mysen [2^ . We note also 
that K c i remains almost constant over the entire concen- 
tration range of interest suggesting that the conversion 
between species is not favoured in the stressed rigid side 
of the glass formation range. On the other hand, the 



TABLE II: Glass optimum condition / = for different SICA 
steps and different cluster types in the case of random bond- 
ing. 



SICA step I 


1 


2 


3 


4 


x c 


0.500 


0.739 


0.690 


0.715 


x c (dendritic) 


0.500 


0.739 


0.667 


0.681 



rapid variation of K c2 can arise from the floppy nature 
of the backbone, when NBO/T > 1.25. The change in 
equilibrium constant following the underlying elastic na- 
ture of the network has been questioned by Eckert and 
co-workers in phosphorus chalcogenides [43| . 



B. Rigidity transitions from cluster probabilities 

From the constraint counting applied on the set of clus- 
ter of a given size I, it is possible to extract a certain 
number of quantities of interest in the context of rigid to 
floppy transitions. 

The simplest case allowed in the case of the SICA con- 
struction is the case of random bonding when no selective 
rules are given on the energies to retain specific clusters 
as starting point for the next cluster generation. This 
means hat the cluster probabilities p n are only given by 
their statistical weights g{E n ) and reduces for instance 
the probability of a CaSiO^ cluster (see equ. (0J) at step 
I = 2 to: 



P2 = PCaSiOs 



x(\ — x) 
(2 - x) 2 



(13) 



or the probability P3R of a three- membered ring SizO§ 
at I = 3 to: 



P3R 



36(1 - xf 



(11 - 10a;)(2 - xf 



(14) 



It is found (see Tabic [HJ that the vanishing of the num- 
ber of floppy modes / occurs at a somewhat higher value 
x c than the macroscopic Maxwell value of 0.50. Also, if 
one considers only dendritic clusters (i.e. no allowance 
for ring structures) , the threshold x c is lowered for steps 
I = 3 and 1 = 4. No change occurs for step I = 2 as 
there are no rings created at this step. The lowering 
of x c for dendritic clusters arises from the fact that the 
stressed rigid ring structures which have a high statis- 
tical weight g(E n ) due to their higher connectivity are 
now absent leading to a number of constraints computed 
from equ. JSJl that is lowered. The probability of finding 
floppy, isostatically rigid and stressed rigid clusters can 
be computed for each SICA step and it is found that for 
dendritic clusters there is a global increase of the floppy 
and isostatically rigid clusters with increasing Ca concen- 
tration whereas the probability of stressed rigid clusters 
is continuously decreasing (see Figure ITU|l . However, no 
threshold is observed as the one suggested by the Ra- 
man analysis and one obtains a single solution / = 
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Ca concentration 



FIG. 10: Probability of finding isostatically rigid and stressed 
rigid clusters as a function of calcium concentration x for 
SICA steps I — 3 (soldi line) and I = 4 (broken line) in the 
case of random bonding. The vertical line corresponds to the 
mean-field rigidity transition I = 4 defined by / = 0. 

corresponding to the glass optimum point defined by the 
vanishing of the number of floppy modes. 



FIG. 11: Probability of finding stressed rigid, isostatically 
rigid and floppy clusters as a function of calcium concentra- 
tion x. The thick solid line corresponds to the basic SICA 
step 1 = 1 and I = 2, the thin solid line to I = 3 and the 
broken line to I — 4. The vertical lines define the intermedi- 
ate phase of width Aa; at setp 1 — 4. The present situation 
for I = 3 and I — 4 corresponds to self-organization with ring 
structures. Notice the kink around x = 0.40 corresponding to 
the stress transition (see text for details). 



C. Boolchand intermediate phase 

We now turn to self-organization that permits to ob- 
tain an elastic intermediate phase |23(,|33(. This elas- 
tic phase has been first observed by Boolchand and co- 
workers in chalcogenides [22 7 but as rigid to floppy tran- 
sitions have been also detected in oxide glasses [44| , there 
is no reason why this intermediate phase should not ex- 
ist in the present system as well. We consider here the 
strutural possibilities that can lead to an intermediate 
phase. It is achieved with SICA by selecting the path- 
ways of cluster production. One starts for example with 
a undcrconstrained (floppy) cluster of size I which exists 
at high calcium concentration. Agglomeration of basic 
units I = 1 onto this cluster is only allowed if the cre- 
ation of a stressed rigid region can be avoided on the 
generated cluster of size (I = 1). This would happen 
if one starts to connect two SiOi/2 tetrahedra together, 
involving an energy gain of E s . On the other hand, if 
x is decreased one will accumulate with this rule iso- 
statically rigid regions on the generated clusters as the 
only allowed connections are either floppy CaO — CaO 
or CaO — SW4/2 bondings. Alternatively, if one starts 
from the low calcium side, self-organization can be ob- 
tained by selecting along the same scheme stressed rigid 
and isostatically rigid connections and excluding system- 
atically the possibility of floppy CaO — CaO bondings in 
the SICA construction. 

With decreasing calcium content, one will be able to 
maintain that rule up to a certain point beyond which 
connections between two silicon tetrahedra (or a stressed 



rigid connection) can not be avoided anymore. The lat- 
ter point corresponds to a stress transition [33| and ap- 
pears only if some medium range order (MRO) made of 
rings is accepted in the construction. The concentration 
range bounded on its low calcium side by the stress tran- 
sition and on its high calcium side by the vanishing of the 
number of floppy modes, defines the intermediate phase 
of width Ax. Results of self-organization are displayed 
in Fig. El The simplest case for self-organized clusters 
is again the case where rings are removed from the con- 
struction corresponding to dendritic clusters which would 
correspond in the limit I — > 00 to Bethe lattice solutions 
|45j or Random Bond Models [l^j obtained in the context 
of floppy to rigid transitions. We obtain here a single 
transition for all SICA steps I either at the mean-field 
value x — x c — 0.50 or closeby (x c = 0.509(7) for 1 = 4). 
No intermediate phase is obtained. The probabilities of 
floppy, isostatically rigid and stressed rigid clusters can 
be displayed as a function of calcium content (Fig. 
and show that the abundance of isostatically rigid clus- 
ters is maximum at the threshold defined by / = 0. This 
is obviously the case for the I = 2 case, but also for larger 
SICA steps. 

The intermediate phase shows up if a certain amount 
of medium range order (MRO) is allowed. This is realized 
in the SICA construction by generating in the construc- 
tion ring clusters such as the one displayed in Fig. 
The requirement of self-organization in the cluster con- 
struction still holds for dendritic stressed rigid structures 
made of at least two connected Si04/2 tetrehadra which 
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Fraction of MRO at x=0.50 

FIG. 12: Width of the intermediate phase at step I = 4 in 
(1 - x)Si0 2 - xCaO gl function of the fraction 

of medium range order (rings) computed at x — 0.50. The 
insert represents the computed number of floppy modes for 
step I — 3 (solid line) and I = 4 (broken line). 

propagate stress in the structure. But now cyclic struc- 
tures such as rings are preserved from self-organization. 
Two transitions are obtained for every SICA steps I > 2 
(see Fig. ^2 The hrst transition one lies always around 
the concentration x C 2 = 0.50 calcium and corresponds to 
a rigidity transition where the number of floppy modes 
vansihes. The second transition that emerges with in- 
creasing MRO is located for I — 4 at x c \ — 0.392 and 
corresponds to the stress transition. When starting from 
a floppy network at high calcium concentration, the pro- 
gressive stiffening of the network can be accomplished by 
requirement of self-organization leading to the accumu- 
lation of isostatically rigid regions and stressed rigid ring 
structures. This will work for any decreasing x up to 
x C 2 below which stressed rigid bondings outside of ring 
structures can not be avoided anymore. 

The two transitions x c \ and x C 2 will define an inter- 
mediate phase Aa; that depends on the fraction of MRO 
allowed and we show that Aa; is anincreasing function of 
the MRO (Fig ll2|) . Furthermore, as the rigidity transi- 
tion / = at x C 2 is almost not affected by the presence 
of MRO, the increase of the width with the latter quan- 
tity arises mostly from the decrease of x c \ with growing 
MRO. Finally, as seen from Fig. ^2 there is kink observ- 
able at x c \ which would produce in first-order derivatives 



such as the energy |47l[ a jump suggesting a first-order 
stress transition at x c \ |23j , |33| and a continuous rigidity 
transition &tx C 2- 

The elastic nature of the network can be also analyzed 
within this framework. Fom Fig. ^] one can see that the 
probability of finding isostatically rigid clusters is max- 
imum in the window Aa;. It is equal to 1 for the I = 2 
SICA step and about 0.5 or 0.6 for the larger steps thus 
providing evidence that the molecular structure of the 
network in the window is almost stress-free. 



VI. SUMMARY AND CONCLUSIONS 

We have shown in the present work from Raman mea- 
surements that a particular transition was observed at 
around 47% calcium in the calcium silicate system, in 
harmony with previous findings [13J . We have elucidated 
the origin of the observed threshold as being the manifes- 
tation of a global softening of the glass structure that sig- 
nals ultimately in a rigid to floppy transition. Constraint 
counting algorithms applied on size increasing clusters 
permit to refine the picture by predicting the possibility 
of an intermediate phase in the Raman threshold region. 

Recent calorimetric and spectroscopic studies have 
shown that this kind of elastic transitions could be found 
in alkaline silicates as well, with a well-defined signature 
of the intermediate phase in both sodium |48| and potas- 
sium glasses ^ij . We are confident that a similar generic 
behaviour should be expected in the alkaline earth sili- 
cate glasses. A definite probe of the existence of the in- 
termediate phase in the present system could be obtained 
using MDSC measurements during the glass transition. 
The results present always a vanishing of the kinetic de- 
pendent heat flow inside that phase |22| and provide an 
unambigous signature for both transitions, rigidity and 
stress. Unfortunately, the glass transition temperatures 
of calcium silicates are too high to be accepted from the 
actual performances of the MDSC apparatus [50j. 
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